Systematic uncertainties of hadron parameters obtained with QCD sum rules 
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We study the uncertainties of the determination of the ground-state parameters from Shifman- 
Vainshtein-Zakharov (SVZ) sum rules, making use of the harmonic-oscillator potential model as an 
example. In this case, one knows the exact solution for the polarization operator II(p), which allows 
one to obtain both the OPE to any order and the spectrum of states. We start with the OPE for 
n(/i) and analyze the extraction of the square of the ground-state wave function, R oc |\Po(r = 0)| 2 , 
from an SVZ sum rule, setting the mass of the ground state Eq equal to its known value and treating 
the effective continuum threshold as a fit parameter. We show that in a limited "fiducial" range of 
the Borel parameter there exists a solution for the effective threshold which precisely reproduces the 
exact IT(/x) for any value of R within the range 0.7 < R/Ro < 1.15 (Ro is the known exact value). 
Thus, the value of R extracted from the sum rule is determined to a great extent by the contribution 
of the hadron continuum. Our main finding is that in the cases where the hadron continuum is not 
known and is modeled by an effective continuum threshold, the systematic uncertainties of the 
sum-rule procedure cannot be controlled. 

PACS numbers: 11.55.Hx, 12.38.Lg, 03.65.Ge 

1. INTRODUCTION 

A QCD sum-rule calculation of hadron parameters [l[ involves two steps: (i) one calculates the operator product 
expansion (OPE) series for a relevant correlator, and (ii) one extracts the parameters of the ground state by a 
numerical procedure. Each of these steps leads to uncertainties in the final result. 

The first step lies fully within QCD and allows a rigorous treatment of the uncertainties: the correlator in QCD is 
not known precisely (because of uncertainties in quark masses, condensates, a s , radiative corrections, etc), but the 
corresponding errors in the correlator may be systematically controlled (at least in principle). 

The second step lies beyond QCD and is more cumbersome: even if several terms of the OPE for the correlator 
were known precisely, the hadronic parameters might be extracted by a sum rule only within some error, which may 
be treated as a systematic error of the method. It is useful to recall that a successful extraction of the hadronic 
parameters by a sum rule is not guaranteed: as noticed already in the classical papers P, 0] , the method may work 
in some cases and fail in others; moreover, error estimates (in the mathematical sense) for the numbers obtained 
by sum rules may not be easily provided — e.g., according to [l|, any value obtained by varying the parameters in 
the sum-rule stability region has equal probability. However, for many applications of sum rules, especially in flavor 
physics, one needs rigorous error estimates of the theoretical results for comparing theoretical predictions with the 
experimental data. Systematic errors of the sum-rule results are usually estimated by varying the Borel parameter 
and the continuum threshold within some ranges and are believed to be under control. 

The goal of this paper is to study systematic uncertainties of the sum-rule procedure in detail. To this end, a 
quantum- mechanical harmonic-oscillator (HO) potential model is a perfect tool (see also Q): in this model both the 
spectrum of bound states (masses and wave functions) and the exact correlator (and hence its OPE to any order) are 
known precisely. Therefore one may apply the sum-rule machinery for extracting parameters of the ground state and 
check the accuracy of the extracted values by comparing with the exact known results. In this way the accuracy of 
the method can be probed. 

We show that the knowledge of the correlator in the limited range of the Borel parameter is not sufficient for a 
reliable extraction of the ground-state characteristics from the sum rule, even if the mass of the ground state is known. 
One should also know the continuum contribution to the correlator with a good accuracy. 

In connection with this observation, we indicate two dangerous points in a typical sum-rule analysis: 

(i) A simple modeling of the hadron continuum by a constant effective continuum threshold leads to uncontrolled 
errors in the extracted hadron parameters. This occurs even in the case when the true effective continuum threshold 
may be well approximated by a constant, as it happens in the HO model considered. 

(ii) The independence of the extracted ground-state parameter of the Borel mass does not guarantee the extraction 
of its true value. 



2. THE MODEL 

We consider a non-relativistic potential model with the HO potential 

V(r) = ^-, r=\r\, (2.1) 
and study the polarization operator II (E) defined by 

U(E) = (27r/m) 3/2 (f f = Q\G(E)\n = 0), (2.2) 
with G(E) the full Green function of the model, 

G{E) = (H-E)-\ H = H + V(r), H Q =f/2m. (2.3) 
The full Green function satisfies the Lippmann-Schwinger operator equation 

G^ 1 (E) = G^ 1 (E) + V, with Go(E) = (H -E)-\ (2.4) 
which may be solved perturbatively: 

G(E) = G (E) - G (E)VG (E) + G (E)VG (E)VG (E) + ■■■ . (2.5) 
For the polarization operator given by a dispersion representation 

K(E) = J J^PM, (2-6) 

the Borel transform [l| has the form 

n(/i) = J dzexp(-z/p)p(z). (2.7) 

Therefore the Borel transform corresponds to the evolution operator in the imaginary time 

HO*) = (27r/m) 3/2 (f f = 0| exp(-H/n)\fi = 0). (2.8) 
For the HO potential (|2.1[) . the exact n(/i) is known [2J: 

^-(sinl^) 372 - (2J) 
Expanding this expression in inverse powers of p, we get the OPE series for n(/i): 
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(2.10) 



and higher coefficients may be obtained from (12. 9|) . Each term of this expansion may be also calculated from 
and (|2.5p , with Ho corresponding to Go : 

oo 

n o(^) = J dzpo(z) exp(-z/p), p (z) = -^=V^ ( 2 - n ) 
o 

The "phenomenological" representation for n(/i) is obtained by using the basis of hadron eigenstates of the model, 
namely 

oo 

Il(n) = J2 Rn exp(-£„/M), (2-12) 

n=0 

with E n the energy of the n-th bound state and R n given by 

Rn = (27r/m) 3 / 2 |*„(r = 0)| 2 . (2.13) 
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The quantity R n determines the square of the leptonic decay constant of the n-th bound state. 
For the lowest states one has 1 

Eo = ^u, E l = 1 ~u, .... (2.14) 

and 

Rq = 2\/2cj 3/2 , R x = 3\/2w 3/2 , .... (2.15) 
For later use we isolate the contribution of the ground state and write 

oo 

U(p) = R exp(-E /p) + U cont {p), n cont (^) ee J dz p phcn (z)exp(-z/p), (2-16) 

Zcont 

where n cont describes the contribution of the excited states (the model has purely discrete spectrum, but we use the 
QCD terminology and refer to the excited states as the "continuum"), z cont is the continuum threshold and p p hcn(z) is 
the spectral density corresponding to excited states. For the HO potential, the continuum threshold lies at z cont = ^u>. 



SUM RULE 



The sum rule claims the equality of the correlator calculated in the "quark" basis (|2. 10|) and in the hadron basis 
(EH: 



i?oe~ BoAl + / dzp phcn (z)e~ z ^ = dzp (z)e- z /» + v 3 / 2 
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(3.1) 



Following we use explicit expressions for the power corrections, but for the zero-order free-particle term we use 
its expression in terms of the spectral integral (|2.1ip . The reason for this will become clear in few lines. 

Let us introduce the effective continuum threshold z c g(p), different from the physical /i-independent continuum 
threshold z cont , by the relation 



n c0 nt(^) = J dzp phcn (z) exp(-z/p,) = I dzpo(z) eacp(-z/fi). 

Zcont Zetl(n) 



(3.2) 



Generally speaking, the spectral densities p p hen(z) and po(z) are different functions, so the two sides of (|3.2[) may be 
equal to each other only if the effective continuum threshold depends on p. In our model, we can calculate n cont 
precisely, as the difference between the known exact correlator and the known ground-state contribution, and therefore 
we can obtain the function z c g(p) by solving (|3 . 2[) numerically. In the general case of a sum- rule analysis, the effective 
continuum threshold is not known precisely and is one of the essential fitting parameters. 
Making use of l|3.2p . we rewrite now the sum rule (|3.ip in the form 



R exp(~E a /p) = U(p,z oS (p)), 
where the cut correlator II(/i, z e s (/-«)) reads 



Il(p, Zcff (/•*)) = -7= / dzi/zexp(-z/p) 
V 71 " J 



3/2 
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(3.3) 



(3.4) 



As is obvious from (|3.3[) . the cut correlator satisfies the equation 

d 



d(l/n) 



logTI(/X, Zes(p)) = Eq. 



(3.5) 



1 Note that, due to the non-relativistic nature of our HO model, the states corresponding to orbital excitations do not contribute to (12.21 
and therefore the excited states contributing to (12.91 are separated in energy by multiples of 2u) from the ground state. 
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The cut correlator is the actual quantity which governs the extraction of the ground-state parameters. It might be 
useful to notice that the relative weight of power corrections in the cut correlator II(/i, z e g(fi)) has been considerably 
increased compared to the initial Hope(h)'- whereas in Hope(^) power corrections are suppressed as l/p 2 compared 
to the free-particle term, in z e s(p)) they are suppressed only as l/y 7 /! compared to the cut free-particle term. In 
the problem under discussion this makes no difference since the power corrections are known precisely. In QCD this, 
however, leads to additional uncertainties since condensates are not always known with good accuracy. 

The sum rule (|3.3[) allows us to restrict the structure of the effective continuum threshold z c g(/x). Let us expand 
both sides of (|3.3[) near tu/fi — 0. The l.h.s. contains only integer powers of u>/fJ,, i.e., even powers of y/iJJJi. Power 
corrections on the r.h.s., on the contrary, contain only odd powers of \J lo j [i. In order that both sides match each 
other, the effective continuum threshold should be a power series of the parameter \Jioj 1 \x: 



LO LO 

Zq + zni— + z 2 — 



(3.6) 



Inserting this series in (|3.3|) and expanding the integral on its r.h.s., we obtain an infinite chain of equations emerging 
at different orders of yjuj/fj,. The equations corresponding to the odd powers of -iJuj/fx do not contain the parameters 
Eq and i?o and constrain the odd-number coefficients z 2 i+\ which provide the cancellation of power corrections on 
the r.h.s. of ()3.3|) . The three lowest-order equations read 

bJZQ 

Ro = ^=z 3/ V/ 2 = J Po (z)dz, (3.7) 
o 

a = A (3-8) 



UJZn 

RqE = -^z 5/ V/ 2 - - _ (z 2 + 4z z 2 ) = / p (z)zdz- - _ y= (zl + iz z 2 ). (3.9) 



OV7T 2^TT^ZQ J 2^TTy/Z 



What is essential is that the i-th equation contains only the variables Zq, . . . , Z\. 

Setting Eq — ^lo and Rq — 2v / 2^ 3 ^ 2 , the equations above yield the following solution for the exact effective 
continuum threshold in the HO model: zq = 2.418, Z\ = 0.142, z 2 = —0.081, etc. 

The following comments are in order here: 

1. Equation (|3.8p rules out the fj,- independent solution z c s = const. 



2. For Eq and Rq within a broad range of values < Rq < i? upP er there exists a solution z c g(fj,, R , E ) which 
exactly solves the sum rule (|3.1[) . Here, the upper boundary i? uppcr is determined from the condition that the 
ground state fully saturates the correlator at /z m i n , the lower boundary of the considered /i- interval: 

-^u PP cr 

exp(— E //j min ) = nnPF,(/imiTi). (3.10) 

For /i min -> 0, i? U pper 2V2w 3 / 2 . 

Therefore, in a limited range of /i the OPE alone cannot say much about the ground-state parameters. What 
really matters is the continuum contribution, or, equivalently, z e g(fi). Without constraints on the effective 
continuum threshold the results obtained from the OPE are not restrictive. 2 

3. The approximate extraction of Eq and Rq worked out in a limited range of values of fi becomes possible only by 
constraining z e g(fj,). If the constraints are realistic and turn out to reproduce with a reasonable accuracy the 
exact z c ff(fi), then the approximate procedure works well. If a good approximation is not found, the approximate 
procedure fails to reproduce the true value. Anyway, the accuracy of the extracted value is difficult to be kept 
under control. 



2 The expected sensitivity of the method should not be overestimated: Imagine, e.g., that we modify the potential as follows: V(r) — > 
V (r) exp(— r/ro). Then the discrete spectrum of states is replaced by a continuous spectrum. However, for sufficiently large values 
TO 3> the power corrections remain numerically almost unchanged. So the Borel-transformed OPE is not very sensitive to the 

dynamics at long distances. 
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Fig. 1: (a) The exact polarization operator (b) The accuracy of the OPE: the quantity (II(/i) — IIope(a 1 )) /IT(/i), where 

n(/i) is the exact correlator and IIope is the result of the OPE involving the first three power corrections. 



The last conclusion is quite different from the results of QCD sum rules presented in the literature (see e.g. the review 
[H). In the next section we shall demonstrate that a typical sum-rule analysis contains additional explicit or implicit 
assumptions and criteria for extracting the parameters of the ground state. Whereas these assumptions may lead to 
reasonable central values, the accuracy of the extracted parameters cannot be controlled. 



4. NUMERICAL ANALYSIS 



In practice, one knows only the first few terms of the OPE, so one must stay in a region of /i bounded from below to 
guarantee that the truncated OPE series reproduces the exact correlator within a controlled accuracy. The "fiducial" 
[l| range of /i is the range where, on the one hand, the OPE reproduces the exact expression better than some given 
accuracy, e.g., within 0.5%, and, on the other hand, the ground state is expected to give a sizable contribution to the 
correlator. If we include the first three power corrections, III, II2, and II3, then the fiducial region lies at ujj [i < 1.2 
(see Fig. [T]). Since we know the ground-state parameters, we fix lu/ fi > 0.7, where the ground state gives more than 
60% of the full correlator. So the working range is 0.7 < uj/^i < 1.2. 

If one knows the continuum contribution with a reasonable accuracy, one can obviously extract the resonance 
parameters from the sum rule (|3.1j) . We shall be interested, however, in the situation when the hadron continuum is 
not known, which is a typical situation in heavy-hadron physics and in studying properties of exotic hadrons. Can 
we still extract the ground-state parameters? 

We shall seek the (approximate) solution to the equation 

00 

Rexp(-E/p) + J dzp (z)exp(-z/n) = ±T pe(m) (4-1) 

in the range 0.7 < ut/p < 1.2. Hereafter, we denote by E and R the values of the ground-state parameters as extracted 
from the sum rule (|4.ip . The notations Eq and Rq are reserved for the known exact values. 



A. /i-dependent effective continuum threshold 



As already explained, since the continuum contribution to the correlator is positive, for any R within the range 
< R < i?upper there exists a solution z c ff(/i, E, R), which exactly solves the sum rule (|4.ip . Clearly, for different E 
and R one has a different, specific continuum contribution n cont (/i, E, R). Thus, without measuring n cont or imposing 
constraints on it based on some other considerations we cannot extract the ground-state parameters! 3 



This is a typical situation when one studies the existence of exotic states, like tetra- or pentaquarks, with QCD sum rules: in this case 
the relevant continuum is not known, and from our point of view, the positive or negative answer to the question whether these states 
exist or not depends mainly on the model used for the continuum. 



6 



In some cases the ground-state energy may be obtained, e.g., from the experiment. However, fixing the ground-state 
energy E equal to its known value E does not help much: for any R within the range < R < R U p P ci- on e can still 
find a solution z c g(fi, R) which solves the sum rule (|4.ip exactly. 

Let us therefore consider constraints on the effective continuum threshold. It is natural to require z e g(p,) > Eq for 
all \i. Then the sum rule (|4.1[) may be solved for any R within the range 0.7 < R/Rq < 1.15. The solution z e g(n) for 
the boundary values of this interval, and the corresponding E{p) = — j- logII(/i, z c g(fi)), R(n), and n cont given by 
(13. 2[) are shown in Fig. [2] Clearly, n cont corresponding to different values of R differ very strongly. 

Fig. [5] also presents the exact effective continuum threshold z c ff(/i) obtained as a numerical solution of the sum rule 
(14. 1|) with the known Eq = |o; and i?o = 2-\/2w 3 ^ 2 , and the corresponding E(/i), R(n), and n cont (/i). 

In the model under discussion, one may expect the exact effective continuum threshold to be somewhere between 
E Q and E^. it is indeed not far from \(Eq + E\) = |w, see Fig. [^a). Requiring, e.g., z e g > \{E + Ex) gives 
0.95 < R/Ro < 1.15, which is also not too restrictive. 



B. Constant effective continuum threshold 



Strictly speaking, the constant effective continuum threshold z c ff(/i) — z c — const is incompatible with the sum 
rule, as it can be seen from Eq. (|3.8[) . Nevertheless, this Ansatz may work well, especially in our model: as can be 
seen from Fig. Ufa), the exact z c ff(/i) is almost flat in the fiducial interval. Therefore, the HO model represents a very 
favorable situation for applying the QCD sum-rule machinery. 
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Fig. 2: Different choices of the effective continuum threshold z c ff(/i) (a) and the corresponding E(n) = — d/dfilog z c s (p)) 
(b), R(fJ,) obtained from the sum rule (14. ip (c), and n cont (/i) given by Eq. (|3.2p (d): 1 [solid (red) line] the exact effective 
continuum threshold as obtained by a numerical solution of (|3.3|l , 2 [long-dashed (blue) line] the effective continuum threshold 
obtained by solving the sum rule (|4.ip for R — 0.7i?o and E = Eq, 3 [dash-dotted (green) line] same as line 2, but for R — 1.15i?o 
and E — Eq, 4 [short-dashed (dark-blue) line] the constant effective continuum threshold z c determined according to Sec. 14 Bl 
In Plot (b), the lines 1, 2, and 3 lie one on top of the other and cannot be distinguished. 
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E{p ,z c ) =E Q , — E(ix,z c )\^ =0, (4.3) 



Now, one needs to impose a criterion for fixing z c . One of the widely used ways is the following [f|: one calculates 

--^-logn^zjE^z,). (4.2) 

The r.h.s. depends on p due to approximating z c g(p) with a constant. Then, one determines po and z c as the solution 
to the system of equations 

d_ 

dp' 

yielding the values z c = 2.454 uj and po/uj — 1, see Fig.[^a,b). The central value of the sum-rule estimate R is obtained 
by setting p — > p , and z c g(p) — > z c in (|3.3p . For this value of z c one has a very good description of II (p) (less than 1% 
deviation in the full range 0.7 < p/u < 1.2) and the stability of E(p 7 z c ) against p is also very satisfactory. Finally, 
in the model under discussion one obtains also a rather good estimate R/ Rq = 0.96, with the function R(p, z c ) being 
extremely stable in the region 0.7 < ui/p < 1.2. Note, however, a dangerous point: the description of II(/i) with 
better than 1% accuracy and the deviation of the E(p, z c ) from Eq at the level of only 1% in the fiducial range leads 
to a 4% error in the extracted value of Rl 

The crucial conclusion from this observation is the following: even when the effective continuum threshold z e e(p) is 
almost flat in the fiducial interval of p, as in our simple model, one still cannot control the accuracy of the extracted 
value of R. As is obvious from Fig.[2jc), it would be incorrect to estimate the error, e.g., from the range covered by 
R when varying the Borel parameter p within the fiducial interval. 



C. Local-duality limit p — > oo 

r— 1 

Let us consider another scheme: a local-duality (LD) sum rule proposed in [7j. This scheme corresponds to the 
limit p — > oo in (|3.3|) and has several attractive features [8J. In the limit /x — > oo all power corrections in the OPE 
vanish and we end up with the simple relation (cf. Eq. I|3.7[) ) 

ZLD 

i? LD = J dzp (z) = gip«& a . (4.4) 
o 

Let us consider the average energy calculated with the cut correlator 

2LD 

dzzpo(z) 



£ld = \- D = |*ld. (4.5) 

dzp (z) 







It is natural to require -Eld = Eq- Then zld = f w and (|4.4|) leads to 



R L u/Ro = P^- 1-05. (4.6) 



As follows from (|3.7p and (|3.9p . the exact values i?o and Eq satisfy the equations 



dzzpo(z) 



n "* / - \3/2 II 3 Z x +4Z0Z2 

Ro= 3^ (a;zo) • So = ^ r 4 ■ (4J) 

dzpn(z) 

Comparing these equations with (14.4)) and (|4.5p . we see that if Eld = Eq, then _Z?ld -Ro- This leads to the 5% 
discrepancy in (|4.6p . Anyway, the estimate (|4.6p is quite good (due to the specific values of the constants z,; in the 
HO model) but its accuracy cannot be controlled. 

Closing this section, we note that the issue of the uncertainties within QCD sum rules (see also [9] for the case of 
light-cone QCD sum rules) appears to be qualitatively similar to what happens in other phenomenological approaches, 
like the constituent quark model, which have indeed many common features with sum rules, as discussed in (To| . 



8 



5. CONCLUSIONS 

We studied the extraction of the ground-state parameters from the polarization operator using various versions 
of sum rules in the case of the non-relativistic harmonic-oscillator potential model. The advantage of such a simple 
model is that both the OPE for the polarization operator and the exact spectrum are known, therefore allowing us 
to compare the results obtained by sum rules with the exact values and to probe in this way the uncertainties of the 
method. 

Our conclusions are as follows: 

• The knowledge of the correlator in a limited range of the Borel parameter fj, is not sufficient for an extraction 
of the ground-state parameters with a controlled accuracy, even if the ground-state mass is known precisely: 
Rather different models for the correlator in the form (|4.1[) — a ground state plus an effective continuum, 
described by an effective continuum threshold z e g (fi) — lead to the same correlator. 

• The procedure of fixing the effective continuum threshold by requiring that the average mass calculated with 
the cut correlator (|3.4|) should reproduce the known value of the ground-state mass [H, @] is, in general, not 
restrictive: a /i-dependent effective continuum threshold z c g(fi) which solves the sum rule (14. ip leads to the cut 
correlator (|3.4[) which automatically (i) reproduces precisely E(pi) = Eo for all values of the Borel parameter /z, 
and (ii) leads to the /i-independent value of R which, however, may be rather far from the true value. 

In the model considered we obtained the following results: 

a. Without constraining z c ff(/i), for any value of R within the range < R/Ro < 1-15 one can find a function 
Zcs(p) which exactly solves the sum rule for 0.7 < to/fi < 1.2. 

b. Requiring 2 eff (/z) > E Q = |w for 0.7 < u/fi < 1.2 gives 0.7 < R/R < 1.15. 

• We studied in detail the standard approximation of the effective continuum threshold with a constant z c . Within 
this approximation, one can tune the value z c by requiring that the average energy E(n) calculated with the 
cut correlator (|3.4|) should reproduce the ground-state energy E in the stability region. In the model under 
discussion, one obtains in this way a good estimate R/Rq = 0.96, with practically /i-independent R. The 
unpleasant feature is that the deviation of R from R turns out to be much larger than the variations of E(n) 
and R over the range 0.7 < u)/ ' ji < 1.2. And, more importantly, error estimates for R cannot be provided. 

• Therefore, we conclude that a sum-rule extraction of the ground-state parameters without knowing the hadron 
continuum suffers from uncontrolled systematic uncertainties (not to be confused with the uncertainties related 
to errors in quark masses, a s , renormalization point, condensates, etc; the latter errors are usually properly 
taken into account). Unfortunately, a typical sum- rule analysis of heavy- meson observables belongs to this class 
of problems: in this case, the hadron continuum is usually not known and is modeled by an effective continuum 
threshold treated as a fit parameter. Then, no estimates of systematic errors for the ground-state parameters 
obtained with sum rules can be given, although the central values may be rather close to the true values. Let 
us also emphasize an important point: as we have demonstrated, the independence of the extracted hadron 
parameters from the Borel mass does not guarantee the extraction of their true values. 

We have nevertheless seen that in the model under consideration the sum rules give good estimates for the parameter 
i?o- This seems to be due to the following specific features of the model: (i) a large gap between the ground state and 
the first excitation that contributes to the sum rule; (ii) an almost constant exact effective continuum threshold in a 
wide range of /i. Whether or not the same good accuracy may be achieved in QCD, where the features mentioned 
above are absent, is not obvious: within the standard procedures adopted in QCD sum rules it is practically impossible 
to control the systematic uncertainties of the obtained hadron parameters. This shortcoming — the impossibility to 
control the systematic errors — remains the weak feature of the method of sum rules and an obstacle for using the 
results from QCD sum rules for precision physics, such as electroweak physics. 
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